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We summarize recent progress in lattice studies of four-dimensional J\f = 4 supersym¬ 
metric Yang-Mills theory and present preliminary results from ongoing investigations. 
Our work is based on a construction that exactly preserves a single supersymmetry at 
non-zero lattice spacing, and we review a new procedure to regulate flat directions by 
modifying the moduli equations in a manner compatible with this supersymmetry. This 
procedure defines an improved lattice action that we have begun to use in numerical 
calculations. We discuss some highlights of these investigations, including the static po¬ 
tential and an update on the question of a possible sign problem in the lattice theory. 


In recent years there has been tremendous progress in lattice studies of four- 
dimensional Af = 4 supersymmetric Yang-Mills (SYM) theory, employing a lat¬ 
tice formulation that exactly preserves a single supercharge Q. 1 ° Advances both 
in the scale of computations and in the algorithms employed^ have provided the 
first ab initio numerical results to be confronted with perturbative and holographic 
predictions for quantities such as the static potential.^ At the same time this work 
has also led to improvements of the lattice construction, in particular the develop¬ 
ment of a new procedure to regulate flat directions in a manner compatible with the 
Q supersymmetry.^ Given the central role of AT = 4 SYM in the AdS/CFT corre¬ 
spondence that relates it to quantum gravity, it is important to continue large-scale 
lattice investigations of the theory away from the regime of weak coupling and for 
arbitrary numbers of colors N. The recent progress, although significant, is only 
the beginning of this effort. 

In this proceedings we briefly review the latest developments and present some 
new preliminary results from ongoing lattice Af = 4 SYM studies. We begin in the 
next section by summarizing the new procedure to regulate flat directions without 
breaking the exact Q supersymmetry. Numerical calculations using the resulting 
improved lattice action exhibit dramatically reduced violations of supersymmetric 
Ward identities and much more rapid approach to the continuum limit. We are now 
carrying out large-scale studies using this improved action. In Section [2] we revisit 
the static potential, checking that the improved action reproduces previous results. 
The improved action has also allowed us to gain new insight into the possible sign 
problem of lattice Af = 4 SYM, which we present in Section [3] We conclude with 
some discussion of other work currently underway and some of the next steps in 
our wide-ranging investigations. 
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1. Supersymmetric deformation of the moduli space 


A convenient starting point for lattice Af = 4 SYM is the direct discretization of 
the continuum action produced by topological twisting, 10 11 
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with repeated indices summed. All indices run from 1 through 5, requiring that the 
theory be formulated on the A\ lattice with five linearly dependent basis vectors 
symmetrically spanning four space-time dimensions.- 2 8 The five complexified gauge 
links U a contain both the gauge and scalar fields, and result in an enlarged V(N) = 
SU(iV) ® U(l) gauge invariance, where N is the number of colors. 

Both the S\J(N) and U(l) sectors necessarily possess flat directions that desta¬ 
bilize the vacuum of the lattice theory. We recently introduced a general method 
that can be applied to regulate these flat directions in a manner compatible with 
the Q supersymmetry. 9 The procedure is to deform the Q {rfD^a term in the 

action, Eq. |T| as 


Q Tr 


V (n) 


£(n 


Q Tr 


V (n) (' 


V[ ) U a (n)+GO{n 


( 2 ) 


where G is a new tunable coupling and 0(n) is an appropriate gauge-invariant 
bosonic operator. This Q-exact deformation modifies the auxiliary field equations 
of motion, d{n) = V^ a J U a (n) +GO(n) Ijv, which determine the moduli space of the 
system. Since Q ij = d, this also modifies the Q Ward identity 


,Tr [Q V (n)]^ = NG (^0( n )j = 0, 


( 3 ) 


where the Tr n)] term vanishes due to the sum over all lattice sites and 

the structure of the finite-difference operator 2? a . 3 ® 

Because the flat directions associated with the U(l) sector produce especially 
severe lattice artifacts,® 9 we take 


0(n) = (det V a b(ri) - 1) = 2Re 'V' (det V ab {n) - 1), 
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where V ab = V£ a is the oriented plaquette in the a-b plane. The corresponding Ward 
identity is ( det ' p ab(n) - 1) =0, implying (Re detV ab ) = 1. However, 

(Im detVab) is still unconstrained, and the SU(iV) flat directions still must be 
regulated as well. We address both of these issues by adding to the action a soft Q 
breaking scalar potential, 
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Previous numerical studies used both this scalar potential as well as an additional 
soft Q breaking potential term that constrained the plaquette determinant. EMI pig. [l] 
illustrates the improvement that results from our new method of lifting the U(l) 
flat directions. Both plots in this figure quantify lattice artifacts by considering 
violations of Q Ward identities, which must vanish in the a/L ^ 0 continuum 
limit. The left plot addresses the Ward identity that fixes the value of the bosonic 
action per lattice site, sb = 97V 2 /2. The right plot shows 1—(Re defP) as discussed 
above; while only the improved action constrains this through a Ward identity it 
must still vanish for the unimproved action as the U(l) sector decouples in the 
continuum limit. In both cases the improved action produces much smaller lattice 
artifacts and approaches the continuum limit much more rapidly, roughly oc (a/L) 2 . 
This scaling suggests that the improved action preserves the Q supersymmetry well 
enough to eliminate any significant effects of dimension-5 operators, all of which are 
forbidden by Q.® 

2. Revisiting the static potential 

The main physical quantity considered by previous studies using the unimproved 
action was the static potential. In our earlier work- -®we found the static potential 
to be coulombic at both weak and strong coupling, V(r) = A + C/r with vanishing 
string tension. For gauge group U(2) our results for the Coulomb coefficient C(X) 
agreed with perturbation theory^ ES f or all A < 2 that we could investigate, where 
A = Ai a t/-\/5 is the continuum-normalized bare ’t Hooft coupling. For N = 3, 
however, we observed a significant departure from perturbation theory for A > 1, 
which could conceivably be the first sign that our results were approaching the 
famous large- N prediction C oc y/X at strong coupling 1 A <C IV. |15 I 16 I 

We are currently repeating and extending these studies using the improved ac¬ 
tion that will reduce any effects of lattice artifacts. As shown in Fig. [2j our pre¬ 
liminary results from this ongoing work are consistent with the earlier observations 




Fig. 1. Continuum extrapolations of Q Ward identity violations on log-log axes with power-law 
fits, comparing the new improved lattice action^ with the unimproved action used in previous 
works .HHU Both deviations of the bosonic action from its exact supersymmetric value (left) and 
(1 — Re detV) (right) are much smaller for the improved action and approach the a/L —»■ 0 
continuum limit much more rapidly, roughly oc (a/L) 2 . For both actions we fix Aj at = 1 and the 
plaquette determinant coupling (G or k), while scaling [i oc 1/L. 
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summarized above. Since the improved action allows us to consider significantly 
stronger couplings, we are in the process of exploring larger A, as well as investigat¬ 
ing N = 4. Comparing results for all three gauge groups U(2), U(3) and U(4) will 
be extremely useful to identify systematic trends, but is numerically challenging® 
due to the rapid growth of computational costs oc N 5 . 


3. New insight into the pfaffian phase 

The improved action allows access to two new sources of information on the possible 
sign problem of lattice AC = 4 SYM. A sign problem is possible because gaussian 
integration over the fermion fields in the lattice path integral produces a pfaffian that 
is potentially complex, J [d4/] T> ' S! oc piT> = |pfP|e*“ where V is the fermion 

operator. All our numerical calculations “quench” the phase e la —> 1 in order to 
take advantage of efficient importance sampling algorithms. ^ For the unimproved 
action subsequent measurements of the pfaffian displayed no indication of a sign 
problem: the phase e lQ exhibited very small and volume-independent fluctuations 
around unity.® 

However, most of our unimproved investigations of the pfaffian phase focused 
on relatively weak couplings, and for volumes larger than 3 3 x 4 we considered 
only Ai a t = 1. As mentioned in the previous section, the improved action allows 
us to study significantly stronger couplings while still controlling lattice artifacts, 
and we are taking advantage of this to begin measuring ( e la ) for larger Ai at < 8. 
Preliminary results from 4 4 lattices are shown in Fig. [ 3 ], and reveal that fluctuations 
in the pfaffian phase increase significantly at stronger couplings Aj a t > 4, apparently 
indicating a sign problem. Part of this increase may result from the larger values 
of the auxiliary coupling G (Eq. [2| that have to be used as Ai a t increases;® larger 
G appears to be correlated with larger pfaffian phase fluctuations at fixed Ai a t- 

The other new development is related to the fact that the Q-invariant deforma¬ 
tion in Eq. [2] breaks an r] —> 77 + clw shift symmetry and lifts the corresponding 
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Fig. 2. Static potential Coulomb coefficients C vs. the bare’t Hooft coupling Ai a t/\/5, comparing 
prelimi nary results obtained using the improved action with previously published unimproved 
results.EEl For gauge group U(2) (left) we observe agreement with perturbation theory. The 
N = 3 results (right) fall significantly below the NNLO perturbative prediction fdQZ] for A > 1. 
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U(l) fermion zero mode. 9 ! The improved action therefore allows us to compute the 
pfaffian of V with fully periodic boundary conditions (BCs). This is not possible for 
the unimproved lattice action, which relies on the fermions’ anti-periodic (thermal) 
temporal BCs to lift this zero mode. In Fig. [3] we compare pfaffian phase measure¬ 
ments with either periodic or anti-periodic temporal BCs for the fermion fields and 
everything else fixed. The contrast is dramatic: anti-periodic BCs produce e la ~ 1 
while fully periodic BCs lead to uncontrolled fluctuations and (e la ) consistent with 
zero. It is not yet clear to us why the pfaffian is so sensitive to these BCs. Even more 
mysteriously, all other observables change very little between these two ensembles, 
despite the apparent presence of a sign problem in one case. 


4. Work in progress and next steps 

Significant work remains necessary to complete both of the ongoing studies discussed 
above. We are still accumulating statistics for several points in each of Figs. [2] and 
[3j For the static potential we are in the process of studying stronger couplings 
and gauge group U(4) to confirm the departure from perturbation theory and ex¬ 
plore whether this behavior is consistent with holographic expectations!^^ We 
are also improving our static potential data analysis, to account for short-distance 
discretization effects when fitting V (r) and to compute renormalized couplings for 
more consistent comparisons with continuum predictions. 
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Fig. 3. Dependence of the pfaffian phase fluctuations on the bare ’t Hooft coupling Ai a t, for 4 4 
lattices and gauge group U(2). 1 — (cos a) vanishes when the pfaffian is real and positive, and its 
large value at strong coupling indicates significant fluctuations in the phase, i.e., an apparent sign 
problem. For each Ai a t we use the smallest acceptable [i and G (Eqs. |5|and|2|) , both of which must 
increase with the coupling. Larger Ai at and larger G both appear to produce larger fluctuations 
in e 101 , potentially compounding the sign problem at strong coupling. 
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For the pfaffian phase, our new results at stronger’t Hooft couplings need to be 
supplemented by measurements on larger volumes V, to see whether the volume- 
independence observed for the unimproved action persists with the improved action. 
We also want to investigate e la for U (N) gauge groups with larger N = 3 and 
4. Both of these goals are challenging;^ the computational costs of our pfaffian 
measurements increase oc N 6 and oc V 3 ~ L 12 . At the same time we continue 
to search for improved qualitative understanding of the system’s insensitivity to 
apparent sign problems at strong coupling or with fully periodic BCs. 

Of course there are many other interesting aspects of lattice Af = 4 SYM that 
we lack the space to discuss in this proceedings. One example is the restoration of 
the other fifteen supersymmetries Q a and Q a t in the continuum limit. This can be 
investigated by measuring violations of discrete R symmetries,® as we previously 
explored with the unimproved action. We are currently revisiting this study using 
the improved action, and will also address the related tuning of a single marginal 
operator that may be required to obtain the correct continuum limit.® ® 

Our highest priority, however, is to compute the anomalous dimension of the 
Konishi operator Ok = XT/ Tr [<f> 7 $ 7 ], where <!> 7 are the six scalar fields of J\T = 4 
SYM. Lattice predictions for this quantity will be complementary to those obtain¬ 
able in perturbation theory,® from holography GSECD or via the conformal bootstrap 
program. 22 This work is advancing rapidly, exploiting both finite-size scaling and 
Monte Carlo renormalization group techniques,® and we hope to report initial 
results in the near future. 
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Fig. 4. Scatter plot of pfaffian phase measurements from independent 3 3 X 4 ensembles using 
either periodic or anti-periodic temporal BCs for the fermion fields (and periodic BCs for all other 
fields and all other directions), with all couplings fixed. While anti-periodic BCs produce e lOL ~ 1, 
periodic BCs lead to uncontrolled fluctuations and ( e XOL ) consistent with zero. 
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